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The effects of fuel burnup 9 fission product poisoning^ and hydrogen 
moderator density variation on reactivity in the AGN-201 reactor are 
considered. A modified one-group perturbation theory is developed and 
applied to changes in parameters resulting from a change in hydrogen 
moderator density* An equation for the reactivity change is obtained 
for three models: a bare cylindrical core 9 a bare core using extrapo- 

lated dimensions 9 and a reflected core. These three equations are then 
used to predict values of the reactivity increase resulting from inter- 
changing a new 3/4 inch thick fuel disk with comparable fuel disks 
presently in the core* 

The results obtained by a digital computer solution of the re- 
activity equations reveal that the increase in reactivity varies from 
0*4392 to 0*7707 ?o 9 depending upon the core model and position of the old 
disk within the core* Because the license of the Naval Postgraduate 
School does not permit a value of excess reactivity above 0*40%* it is 
concluded that a simple interchange of disks in this manner would pro- 
duce too large a value of excess reactivity* 
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NOMENCLATURE 



ENGLISH 


SYMBOLS 


A P 


= molecular weight of polyethylene 


A U 


= atomic weight of uranium 


B 


= buckling of the reactor system 


d 


= extrapolation distance 


D 


= diffusion coefficient of core 


D 

r 


= diffusion coefficient of reflector 


f 


= thermal utilization 


g 


= non-l/v factor 


h 


= distance of nearest edge of fuel disk above or below the zero 
ordinate of the z-axis of the cylindrical fuel core 


H 


- height of fuel core 


H' 


= extrapolated height of fuel core 


I 


= effective resonance integral 



= multiplication factor for an infinite reactor system 
L = diffusion length in core 

= diffusion length in reflector 
M = perturbation operator 
M+ = adjoint of M 



M 9 


= perturbation operator of perturbed reactor system 


M f 


= mass of fuel per disk 


M 

P 


- mass of polyethylene 


N 


= atom density 


N 

a 


- Avogadro e s number 


N 

o 


= atom density at time equal to zero 



s 



p 


= resonance escape probability 


Pm 


= promethium concentration 


PWR - 


reactor power 


R = 


radius of fuel core 


R f = 


extrapolated radius of fuel core 


S - 


samarium concentration 


Soo = 


equilibrium samarium concentration 


t = 


time 


T = 


reflector thickness 


T = 

eff 


effective neutron temperature 


H 

3 

11 


moderator temperature 


T = 
r 


radial reflector thickness 


T « 
z 


end reflector thickness 


v d = 


disk volume 


z d = 


centerline distance from a disk to the zero ordinate of the 
z-axis of the fuel core 



GREEK SYMBOLS 





- fission-product poison yield of promethium 


% 


= reflector savings 


Sh 


= disk thickness 


\r 


- radial reflector savings 


u 


= end reflector savings 


Ao 


= fraction of total change in core moderator mass occurring 
in a fuel disk 




= fast fission factor 


V 


= net neutrons available after fission 


Kp 


= decay constant of promethium 
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'Ur 

V 

f 

f 

T p 

r 

o 335 



, lit 



<r ft 

<*> 

<r, 

£ f 

fo 

0 

0 

0' 



transport mean free path 

average number of neutrons released per fission 
average logarithmic energy decrement per collision 
reactivity 

density of polyethylene 
density of uranium 
density of U-235 
density of U-238 

microscopic absorption cross-section 
microscopic fission cross-section 
microscopic scattering cross-section 
macroscopic absorption cross-section 
macroscopic fission cross-section 
macroscopic scattering cross -sec tion 
Fermi Age 
neutron flux 
average neutron flux 
peak neutron flux 

neutron flux of perturbed reactor system 
importance function 
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I . INTRODUCTION 



The AGN-201 nuclear reactor is a small, safe, low-power, self- 
contained reactor designed and produced by Aerojet-General Nucleonics 
Corporation of San Ramon, California* It is designed primarily for use 
in education and research not requiring a high neutron flux, but featur- 
ing maximum safety, high analytical sensitivity, and ready accessibility 
of components* The fuel core of the reactor is composed of ten circular 
fuel disks of varying thicknesses, piled on top of one another to form 
a cylindrical shaped core* The fuel disks are made of polyethylene and 
have uranium fuel homogeneously interspersed in them. Appendix A con- 
tains a general physical description of the overall reactor and also a 
detailed discussion of its fuel core* 

The AGN-201 nuclear reactor located at the Naval Postgraduate School, 
Monterey, California, has been in licensed operation since April 29, 1957. 
From that date, until January 31, 1962, the reactor was operated at a 
maximum power level of 100 milliwatts* During that time period no de- 
crease in excess reactivity was noted* Thereafter, the license of the 
Naval Postgraduate School was amended to permit continuous operation up 
to 20 watts and intermittent operation at power levels up to 1000 watts* 

On February 3, 1962, the excess reactivity of the reactor was measured 
and recorded as 0*354%* From that time until February 26, 1969, the 
value of the excess reactivity has steadily decreased to a value of 0.281% 
or a drop of 0*073%. 

Several reasons which would account for the observed loss in re- 
activity have been postulated* The first cause of this loss could be 
due to fuel burnup. Since the reactor is highly enriched with uranium-235, 
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(approximately 20 %) , the decrease in reactivity would be proportional to 
the fraction of U-235 depleted by fission,, A second cause of reactivity 
drop could be the buildup of reactor poisons due to fission of the fuel. 
These fission fragments and their decay products have large neutron absorp- 
tion cross-sections o If they are produced in appreciable amounts they 
could effectively alter the neutron balance in the reactor by affecting 
neutron multiplication. Thirdly, the loss in reactivity of the reactor 
might be attributed to a loss in moderator within the reactor core. This 
problem is unique to this type of reactor. Polyethylene, which is the 
neutron moderator, when subjected to radiation, tends to decompose. This 
decomposition is accompanied by the liberation of hydrogen gas. The ef- 
fect of this phenomenon on the reactivity is especially pronounced be- 
cause the hydrogen is bled off from the fuel core at varying time inter- 
vals which reduces the density of the moderator within the core. 

In September, 1968, the Naval Postgraduate School received a new 
fuel core disk for the purpose of adding it to the core in some manner 
to attempt to raise the value of the excess reactivity of the reactor. 
However, the limitation imposed upon this addition is that the resulting 
value of the excess reactivity of the reactor must not exceed 0.40%. The 
new disk has a known fuel enrichment and a higher polyethylene density 
than any of the disks presently in the core. The properties of the new 
disk are presented in Appendix A. 

The objective of this thesis is to determine the change in reactivity 
which would result from interchanging the new fuel disk with a disk of 
comparable size currently in the reactor. To accomplish this objective, 
first an analysis is made of the three possible causes of the reactivity 
loss due to reactor operation to determine which have significant effects, 
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and what changes in reactor properties result from these effects* Second, 
using a modified one-group perturbation theory, an equation for the re- 
activity as a function of these changes is derived for three models des- 
cribing the reactor core* These models treat the core separately as a 
bare cylinder, a bare cylinder having extrapolated dimensions, and a re- 
flected cylinder* Third, a digital computer program is used to solve the 
reactivity equations using assumed changes in reactor properties. Since 
the actual value of reactivity drop due to operation is known, the core 
properties at present can be found* Using these results and the proper- 
ties of the new fuel disk, the reactivity equations derived from perturba- 
tion theory are solved by a second computer program to yield values of the 
excess reactivity caused by the perturbation of replacing an old disk with 
the new disk. 
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II a VARIABLES WHICH AFFECT REACTIVITY 



A. FUEL RURNUP 

The first of the three probable causes of decrease in excess reactiv- 
ity to be considered is that of fuel burnup. From published literature 



actor operation (24 hours a day, 7 days per week) at a 100 milliwatt power 
level would result only in a burnup rate of approximately 30 micrograms/ 
year. On the basis of this information it can probably be anticipated 
that even at higher power levels, the contribution of fuel burnup to a 
reactivity drop will be slight, if any. 

Since the reactor is fueled with only a single fissile substance, 
U-235, and there are no fertile materials present, (it is assumed that 
due to the high enrichment of U-235 there is no compensation for fuel 
depletion through the production of fissile Pu-239 from neutron absorp- 
tion by U-238), the depletion rate of the fuel is governed by the equa- 
tion 



tion cross-section, and 0 is the neutron flux. Rearranging and 
integrating (1) yields 



concerning the AGN-201 reactor [1]\ it is predicted that continuous re- 




( 1 ) 



where N is the atom density of the fuel, 0^ is the thermal absorp- 




' £ dt 



The integral exponent 




o 



, can be expressed as t , 



where 0 i s a n average value of the flux over time t. Then 



1 



Numbers in brackets refer to Bibliography. 
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( 2 ) 



From recorded data contained in the operators log kept at the reactor 
facility, as of December 9, 1968, the reactor had been operated approxi- 
mately 8*69 x 10^ watt-min. Although operation occurred at power levels 
varying from 0.1 watt to 1000 watts, 20 watts was most frequent, and 
therefore was assumed to be a representative average value. This gives 
an average time of operation of 

t = 2.61 x 10^ sec. 

An equation given by Lamarsh [5] relating reactor power and neutron flux, 
assuming a recoverable energy of 200 Mev per U-235 fission is 

3> - (PWRXT^Q 1 ^* ip 13 

C7.i9xn f )c^) 



where PWR is the reactor power, i s the effective neutron tempera- 

ture, Mj. is the mass of fuel, and g^ is the non* j /v factor. Since 

-6 

PWR = 20 x 10 megawatts 

T = 340° K 
ef f 

g f (T) = 0,9269 

and = 0.667 Kg 

evaluation of equation (3) yields 

— 8 2 
0 = 8.3 x 10 neutrons/cm -sec 



and the quantity 

« 1.16 x 10‘ 6 (4) 

where 

^ = THa "Vivff 1 ^ 3* 

Substituting (4) into (2), 

7f 0 ’ e - 0.9999988 
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Therefore it is concluded that fuel burnup is negligible at the assumed 
average power level of 20 watts* 

B. FISSION PRODUCT POISONING 

As previously mentioned , reactor poisons, having large neutron absorp- 
tion cross-sections, have an adverse effect on reactivity. Their chief 
influence on reactivity is through the multiplication factor by decreas- 
ing the thermal utilization* In discussing reactor poisons in the AGN- 
201, it is felt that in order to affect reactivity, a poison produced in 
this particular reactor must meet three criteria* (1) It must have a 
high cross-section* (2) It must have a solid form* (3) It must be 
stable. Although as many as 200 species are apt to be produced during 
the fission process, only two have large enough cross-sections and are 
produced in sufficient quantity to warrant attention; namely Xe-135 and 
Sm-149. 

In analyzing these two, it is seen that Xe-135 meets neither condi- 
tion (2) nor (3) above* Being an inert gas, it would be removed from 
the reactor at the time when the hydrogen buildup is bled off from the 
core. Thus the decrease in reactivity would be a function of time, being 
zero at a time equal to zero, buildup to some value, and then return to 
zero at the time of the hydrogen and poison removal* However, this is 
contrary to observed evidence of reactivity behavior. In addition, Xe- 
135 is not stable but is an intermediate isotope in the chain 



Because it is not a stable isotope, any xenon poison produced during 
operation would eventually decay after reactor shutdown, resulting in 




* Bo 



(stable) 
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no long term effect on the reactivity. Hence, Xe-135 can be effectively 
ruled out as a cause of long term reactivity drop, 

Sm-149, on the other hand, meets all three requirements. It has a 
thermal absorption cross-section equal to 58,500 barns at 20°C. It is a 
rare earth metal and is the stable decay product of the chain 
Nd * Pm 1 Om (stable). 

The neodynium decays relatively quickly compared to the promethium so that 
for practical purposes the Pm-149 may be assumed to be produced directly 
as a result of fission. The concentration of Pm-149 is then governed by 
the equation 

* p £,0-A p Pm (5) 

where Pm is the promethium concentration, tfp is the fission-product 
poison yield for promethium, is the fission cross-section of the 

fuel, 0 is the neutron flux, and Ap is the decay constant for prome- 
thium. 

The rate of change of the concentration of Sm-149 can be expressed 
as 

= ApPm -<r o 5 0S (6) 

since samarium is stable and can only be lost through neutron absorption. 
In equation (6) S is the samarium concentration and is the absorp- 

tion cross-section of samarium. Combining (5) and (6) and integrating 
yields an expression for Sm-149 concentration as a function of time. 



Since 



= 0.0113 



£.f = 0.0646 cm 

_ * oc lrt '20 2 

T 0 = 5^85 x 10 cm 

Ap = 3o56 x 10 ^sec ^ 

and from the previous section it was assumed that 

t = 2.61 x 10^ sec 

8 2 

0 = 8.3 x 10 neutrons/cm -sec. 



Then 



S = 1.313 x 10 12 cm" 3 . 

It can be seen from equation (7) that as t — * 00 , S would no longer be 

a function of time but would approach a constant "equilibrium concentra- 
tion". For purposes of comparison let t *°° in (7). Therefore, 



c-s ^ ^ _ 16-3 

= ^TT - 1.178 x 10 cm 

Thus it is seen that the concentration at present is less than the 

4 

equilibrium value by a factor of approximately 10 . The macroscopic 
cross-section of the samarium is 

= ScT a 5 * 7.68 x 10" 8 cm" 1 . 

The change in reactivity due to samarium poisoning can be determined 
from the equation [5] 



f = gq/S* 
0 



where 



V = 2.43 

and p ~ £ « 1. 

Thus , 

$ = -0.48 x 10"So. 

Using the same equations with the equilibrium concentration of Sm-149 
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yields 

j 1 = -0.4397. 

4 

As can be seen, there is a difference of a factor of 10 between the 
reactivity values. Further, using the same equations to solve for the 
time required to reach equilibrium concentration yields a value of 8.11 
x 10 7 sec. This is a factor of 40 times larger than the time of re- 
actor operation to present. 

From this analysis it can be concluded that samarium poisoning does 
exist in the reactor, but it has not reached its equilibrium value. 

Due to the relatively short period of operating time and the low flux level 
produced during operation of the reactor, the decrease in reactivity due 
to samarium is several orders of magnitude smaller than the observed value 
of 0.0737o. Consequently the effect of samarium will be neglected. 

C. MODERATOR CONCENTRATION 

The third possible cause of reactivity drop is that of loss of modera- 
tor. As already noted, the moderating material in the AGN-201 reactor is 
polyethylene. Chemically, the molecular formula of polyethylene is (CH 2 ) n 
indicating that it is composed of long chain (high polymer) molecules 
formed by many (Cl^) groups. 

-ch 2 ch 2 ch 2 ch 2 « 

In fact, the molecular weight of commercial polyethylene is between 
18,000 and 22,000. From a nuclear standpoint, polyethylene is a long 
series of carbon atoms connected in chains with bare protons (hydrogen 
nuclei) attached to them. When high energy neutrons are introduced into 
a polyethylene medium, it is possible for a neutron to collide with the 
proton and break the proton away from the main chain. When this happens, 
cross-linking occurs, whereby new bonds are formed between adjacent high 
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polymer molecules accompanied by the evolution of a hydrogen molecule. 

«ch 2 ck 2 ch 2 ch 2 - -ch 2 ch 2 ch ch 2 - 

AAA/V— » | +H 9 

-CH 2 CH 2 CH 2 CH 2 » -CH 2 CH 2 CH ch 2 * 2 

This is referred to as unsaturation of the polyethylene. The main point 
to be noted is that when polyethylene is irradiated, gas is evolved, the 
majority of which is hydrogen, along with slight traces of methanes. The 
assumption made, upon all that follows in this thesis is that due to 
hydrogen evolution and subsequent removal of the hydrogen from the fuel 
core, the mass of the polyethylene moderator is decreased and that this 
decrease is equal in magnitude to the mass of hydrogen removed. 

From research conducted with the AGN-201 in March ,1962, shortly after 
the reactor license was amended to permit higher power operation, it was 
noted that gas evolution from the core occurs primarily at power levels 
between 200 and 1000 watts. [7] This gas evolution takes place slowly 
and over a period of about two days following operation. No gas evolu- 
tion was observed in relatively short periods of operation at 20 watts. 

A sampling of the gas revealed traces of an inert gas mixed with the 
hydrogen, the decay rate of which corresponded to xenon. The fact that 
the beginning of the observed decrease in excess reactivity coincided with 
the start of the reactor tf s operating at high power levels tends to con- 
firm that hydrogen evolution is the cause of the reactivity drop. 

In the theory that follows, the decrease in density of the poly- 
ethylene moderator manifests itself in two ways. First, it appears 
directly in the moderator density term, $ Secondly, it causes a change 

in the fission and absorption cross-sections, (T ^ and <T of materials 
in the reactor. Because the amount of moderator in the fuel core is 
decreased, the neutrons are not as well moderated and consequently the 
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average energy of the neutrons is at a higher energy value* When this oc- 
curs, the cross-sections of the materials in the reactor, which are 
energy dependent, are no longer those corresponding to 0*025 ev, but rather 
those which correspond to a slightly higher energy® 

The energy at the peak of the energy curve is related to a tempera- 
ture known as the effective neutron temperature® The neutron temperature 
is related to the density of the polyethylene in the following manner. 
Coveyou [9] has established the relation 



where 

and 



T -- = T (1 + 0*46 & ) 
eff m 

A - W 4 *') 

< r. X p Nc, 

' “A 



for moderators with mass from 1 to 25 and ^ from 0 to 1®0® For heavy 
moderators, Cohen [9] gives 



T «£f- T n, <1+0 - 3 *>- 

From a survey of research done in the field of polyethylene, Cooke [3] 
concluded that the effective neutron temperature, T^^, of neutrons in 
polyethylene was between 336 and 347° K® Since it is known for the re- 
actor that A “ 0.2931 and T = 293° K, and T -- was assumed to be 340° K, 

m eff 

the numerical coefficient in the expression for T was modified to a 

eff 

value of 0*5459 to fit the p@lyethylene data® The resulting equation is 
thus , 

T .. - T (1 + 0,5459 L ). (8) 

eff m 

From the equation linking the density of polyethylene to the effective 
neutron temperature, the neutron temperature can in turn be used to find 
the cross-section at that temperature through the following equation, 
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? <r 3 <t) 

Thus, from these equations it can be seen that, the effect of decreasing 
the polyethylene density is to increase the neutron temperature which 
decreases the neutron cross-section of materials in the reactor. 

From work done with the reactor in June, 1963, based on test runs 
of 1000 watt-min, the rate of hydrogen evolution from the fuel core was 
found to be about 0,72 moles/KW-hr, [8] The number of kilowatt-hours 
of operation from February, 1962, to February, 1969, the period during 
which the drop in excess reactivity has occurred, is known to be about 
16,07 KW-hr, Therefore, the number of moles, and consequently the weight 
of hydrogen evolved can be approximated. This results in a value of 
hydrogen weight of 23,14 gm which when equated to a change in polyethy- 
lene mass in the core gives a value which can be used to check against 
that predicted by perturbation theory, 

D a CONCLUDING STATEMENTS 

The conclusions drawn from the preceding sections can be summarized 
as follows. Fuel burnup can be neglected as a cause of reactivity de- 
crease on the grounds that the time of operation of the reactor has not 
been long enough nor is the assumed average level of operation at 20 
watts high enough to cause any significant depletion of the U-235 fuel 
in the reactor core. The two reactor poisons most likely to be produced 
during operation of the AGN-201, xenon and samarium, cannot account for 
the drop in excess reactivity. The Xe-135 is either withdrawn from the 
reactor simultaneously with the hydrogen during bleedoff or decays to a 
more stable isotope during the frequent periods of reactor shutdown. The 
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extremely small concentration of Sm-149 present in the reactor does pro- 
duce a decrease in the excess reactivity but its magnitude is several 
orders smaller than the observed decrease and is negligible,. The de- 
crease in polyethylene moderator due to hydrogen evolution resulting from 
neutron irradiation cannot be ruled out as a cause of reactivity drop* 

It is assumed that this cause is solely responsible for the decrease in 
reactivity* Also, observed evidence seems to indicate that this is the 
most probable cause of the reactivity decrease* The change in polyethylene 
density will appear in the theoretical equations that follow both as a 
direct density change and also as changes in neutron cross-sections of 
reactor materials* This is because the cross-sections are functions of 
neutron temperature which in turn is a function of moderator density. 
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XII* PERTURBATION THEORY 



A. INTRODUCTION 

To analyze the effect on the behavior of a nuclear reactor caused by 
small non-uniform changes such as those discussed in the preceding para- 
graphs , a method known as perturbation theory has been developed. If the 
problem, which has been outlined, involving the AGN-201 were of such a 
nature that uniform changes occurred throughout the reactor, then the re- 
activity could simply be found by recomputing new reactor constants and 
recalculating a new multiplication factor to find the resulting reactivity. 
When non-uniform changes within the reactor occur, as is usually the case 
in practice, then perturbation theory must be employed. Examples of non- 
uniform perturbations include insertion or withdrawl of control rods into 
the reactor, non-uniform fuel burnup, non-uniform fission poison accumu- 
lation or as it is believed in this case, non-uniform removal of modera- 
tor from the core. The assumption upon which perturbation theory is based 
is that the perturbations or changes within the reactor are small and that 
the flux is not distorted substantially in the region of the perturba- 
tion. Since this is assumed to be the case in analyzing the AGN-201, 
perturbation theory provides a convenient method of handling an other- 
wise difficult problem. 

B. MODIFIED CRITICALITY EQUATION 

Perturbation theory is referred to as one-group, two-group, or many- 
group depending upon the type of theory used to describe the energy of the 
neutrons in a critical reactor. The simplest is one-group and will be used 
here. The criticality equation describing a one-group reactor is 
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where 



i * i x B a 

a '’r-F 



Strictly speaking, this equation is used in a one-group perturbation 
analysis. However , it was felt that improved accuracy of the results 
could be obtained by using the Fermi Age-Dif fusion model to include neu- 
trons in the slowing down region. The usual form of the criticality 
equation in this case is 



.0^ o 



1^(6 | 
i t ~ 1 



(9) 



where K 00 *^£p , f , C ~ is the fast non-leakage probability, 

and 0* ^ ^ ) is the slow non-leakage probability term. However, the 
method used to calculate t} 0 , (the Fermi Age equation)^ is not valid in 
the slowing down region for the case of an hydrogeneous moderator be- 
cause the assumption of many collisions by the neutrons does not hold. 

If the moderator contains hydrogen to any degree, it is possible for 
the neutrons to become thermalized in a few or perhaps even a single 
collision. Weinberg and Wigner [9], have developed a single collision 
kernel to account for this type of moderation. The result is that the 
usual fast non-leakage term is replaced by 



* Qn "' (“fe) (10) 



where is the cross-section for this type of interaction and is 



given as 



SL*- 



(3? 0 y* 

where Co is the experimental value of the Fermi Age. Substituting 
(10) into equation (9) gives 

--Jgop.C. lr J. _ | 

I + l/B* “ 1 



T 



This is the form of the modified one-group criticality equation to which 
strict one-group perturbation theory will be applied* 



C* DEVELOPMENT OF PERTURBATION THEORY 



In this section* one-group perturbation theory is employed to derive 
an expression for a change in reactivity from changes in various nuclear 
parameters . 

Let 0 *' B tan V 2 0 J 

The modified one-group criticality equation for an hydrogenous reactor 
developed above becomes 



<«c 



iu % e l - 1 

This equation is the direct result of a modified diffusion equation of 
the following form* 



+ div D grad 0 = O (ii) 

where is the source term* <£$ ft is the absorption term 

and div D grad 0 is the leakage term* Equation (11) can be written as 

where M is the operator 

M = VZPC - £* * d»v D 

Referring to the discussion of the adjoint operator in Appendix B, the 
adjoint equation* which will be of use later in the development but which 
is presented here* is 

fi t 't - o 

M+ is the adjoint of M and the function is referred to as the im- 

portance function* For one-group perturbation theory* M is self-adjoint 
and 0 and P are everywhere proportional* Terms in equation (11) which 
are variables are P, C* j * and D. Considering first perturbations 
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or changes in ? 9 C s 2^ 9 and 2 a 



s> 



P 9 = p + % p 

c 9 = c + $c 

tf ' - £,+ S£ f 

£,’ - 2c+ «« 

The unprimed terms refer to unperturbed quantities , the delta terms are 
magnitudes of the changes or perturbations , and the primed terms refer 
to the resulting perturbed quantities* To return the reactor to critical 
a change must be made in *1) . This is not possible in reality since V 

is a physical constant* However*, the concept of compensation in this man 
ner is useful in the analysis , so that 

V* 5 V + 

For the perturbed reactor 

MV*o 

where 

M' -- (V* + -< + 0 grac* 

Neglecting higher order terms*, since for perturbation theory to be valid 
the perturbations must be small, we get 

tf* s «■ %V€PQ.2 f '“SPveCZf + 

+ tt'f'VzP c- Za, - %to> * oti v 0 <jrac( (12) 

Substituting the quantity M into (12) gives 

M's M + %v*P * ^€CS f + 

This may also be written as 

M' 3 M * A„ 

where 



A 0 • wtPQ.Zjf 



Therefore 



(n* A e )$' s o 

M + ^ s o 



(13) 



(14) 



Multiplying (13) by 0 , (14) by 0 , subtracting (14) from (13) and 
integrating the result over the volume of the fuel core yields 

X KM + Ao)0'dV - X tf'nV olV * o 

Expanding the first integral, 

X S‘M0'ctV ♦ J) VA,0'olV -J' v 0'M , V'dV*O 

Combining the first and third integrals, 

X OM0'- V)dV *i v 0A o 0'dV* O (15) 

However, from Appendix B, by definition of the adjoint operator, 

X ( erics'- 0 " mV) dV'-o 

Therefore, (15) becomes 

S v V-A. 0'dV* O 

Substituting the expression for Ao gives 

X + ^Pz^eC2 f +-*SC^€P2 f + *§2 f V£PC ~12 a )0'dV = O 

Solving for and since S - ~ ^ 

? - X frCSPSCgf + £ F£r I yS±jPC - -y )0 , dV 

X K* p C2f)0'e*V 

Since it is assumed that the perturbation is small;, 

0 

and because the one-group operator is self-adjoint, 

0 « 0 

and the expression for $ becomes 



X C^€PCS f ) 0*dV 



( 16 ) 
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Next, consider a change in D only,* 

D' = D + ‘iD 

Then 

M'-- (v* %i>X^PCZf) + d iv (D+^D^grad 

Following the analysis previously outlined*, the operator for the per- 
turbed system can be written as 

n ' 3 m * A a 

where the perturbation operator is 

A c * * div ID grad. 

This results in a reactivity of 

J _ J, 0ctiv -aO 9rod cW 
J v 0£PC2^)« l dV 

By employing vector identities*, 

S v f div ID grad <0 dV - X div ( 0 ID grad $OdV - X ^0 Cv0VdV 

However, since the flux vanishes at the reactor core surface 

X div (tf'SD grad S»)dV - f A g rad 0 • n dA * O 

Hence , 

X 0 div "5D grad 0 dV = “ J v "SO CV0) z dV 

Therefore, (17) can be expressed as 

S - (18 , 

J y (t^PCSf^cW 

Changes in P*, C s 2^ s 2 a and D occurring simultaneously produce a re~ 
activity represented by the sum of (16) and (18) or 



^ . Jv(l>CaPeCS f + ‘SCe PS t + ^g^PC)-^^ Q l0 Z ’- ^PCV0) 4 }dV 

j v CZf)# 1 - dV 
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D. APPLICATION OF PERTURBATION THEORY 



1 • Bare Cylindrical Geometry 



In this section the expression derived for the change in re- 



activity is solved for the case in which the fuel core is unreflected 
and the actual geometrical dimensions of the core are used. The signi- 
ficant dimensions and coordinate system used are shown in Figure 1. The 
numerator of the expression for the reactivity contains terms which repre- 
sent changes in particular reactor parameters in the perturbed system. 

The perturbation is assumed to occur only in one disk of the core at a 
time and is uniform within that disk. The total reactivity of the core 
is found by summing the reactivities of the ten disks. ">ince the perturba- 
tion is zero everywhere in the core outside the disk under consideration, 
the limits of integration in the numerator are from zero to R in the 
radial direction and from h to h *f ^ h in the vertical direction. The 
denominator, however, contains no perturbation terms and is integrated 
over the entire volume, i.e., from zero to R and from zero to H. 



Let 



I, = c(V 

l z = J v 1>CZJ€P2 f 0^ 

1 3 = J v 

I* = T V 





Therefore , 




( 19 ) 
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z 




Fuel Core Dimensions 
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Evaluation of these I terms is presented in Appendix C. Substituting the 
resulting expressions for the I terms in (19) 9 and cancelling like terms 
results in 

j- - <3 t It p CQ - _}_£ a Q ■ 1>D + (ff )* Y] 

TTx'SPCSf 

Q 

This expression can be further reduced by factoring out the tT term and 
dividing each term in the numerator by the denominator. This step is 
justified only if it is assumed that the values of P 9 C, 2^, and D 

appearing in the numerator* which refer to disk properties 9 do not vary 
significantly from their corresponding values in the denominator which 
refer to the entire fuel core. Making this assumption makes the final 
reactivity equation much more tractable „ Then* 




Since ^ where f 

can be written as 






Vi 



PCS* + ( H ) Q ]^ D } 




and = V 



JGf 



the expression for 



s 



f - 2.(5 1 

5 TT \ P 



3 >c_ vs* 

C 2 c 



I H2q 



I 



<o=,C,£ 






Letting * 

and since * 



e- 



U 



, the final expression for the reactivity for this 



case is 



, tc. l___ Hi* __ L^_e \ 

TT \ P C £ f K W C K„C D ) 
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2. Extrapolated Dimensions 



The derivation of the equation for the reactivity for the case 
using extrapolated dimensions is similar to that of the bare core geo- 
metry. The difference results from the fact that for this analysis the 
flux is no longer considered to be zero at the physical boundary of the 
fuel core but instead goes to zero at a hypothetical distance, d, beyond 
the actual boundary. The resulting flux distribution more accurately 
represents the actual situation within the reactor core since experimental 
evidence indicates that the flux is not zero at the core boundary but has 
a finite value. It can be shown by transport theory that for a planar 
free surface, d has the value 



cl = 0.71 A ir 

The transport mean free path, X^ P is found from the equation 



A 3 



I- ^ 

But since is small for the polyethylene moderator, it is assumed 



that 



K 



tr 






Substituting this into the equation for d gives 



For a curved free surface, such as the outer edge of the cylindrical 
fuel core, the above equation for d can still be used and introduces 
negligible error provided that d is small compared to the radius R. 

This condition will be shown to be true for this reactor. Therefore, 
because of the extrapolation distance, the extrapolated radius and extra- 
polated height of the fuel core become 

R 9 * R -f d 
H 9 - H + 2d 

and the parameters and A3 in the fl rx expression become 
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o< 



3 L .*fQS 

R' 



/6 



T r 



H' 



From this point on^ the development of the reactivity equation follows that 
outlined in the preceding section with two exceptions,, First, several 
of the zero order Bessel functions in the reactivity equation for the 
bare core geometry case were ^ero because of the boundary condition that 
the flux equal zero at R„ When using the extrapolated radius R' these 
Bessel functions are not equal to zero and are retained in the equation. 
Secondly, using the extrapolated height H 0 , rather than H in many of the 
integrals involving trigonometric functions yields a slightly different 
result. Because much of the development of the reactivity equation is 
repetition of the bare core geometry case it will not be repeated here, 
but it can be shown that for this case that 



p _ _L_ / It + ^ L_ •ss* Li b -SO \ 

j " Zt\ v P C £ f <ocC 2o Ko»C D / 



where 



T - Hit! 4 X <; n 

” H' % H' 



Z 






qnh 

H' 



A = 



TTH | 
ZH* 



5 in 



QTTH 

H' 



b = 




(ML' 


K W 






E 


\ H 1 , 


) T 






W = 


F-Sh 


+ ± 


s in 


iZil _ _L 


arr<MSh) 

im — ■— r; 


H' 


Z 


M' Z 


H 



X s JJ (*R) - Jo 0*0 1 * 0*0 



and 



E* T 0 1 OR)* J , : 1 OR) 
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3o 



Reflected Core 



Since the fuel core of the AGN-201 reactor is surrounded by a 
graphite reflector, the expression derived for the reactivity in this 
section probably is the most realistic and accurate of the three cases 
discussed, even if it is only an approximate solution* As it turns out, 
the problem of a cylindrical core with radial and end reflectors, as is 
the case for this reactor, is more than can be handled in an exact analyti 
cal manner* The problem arises from attempting to find a form in which 
to represent the flux distribution so as to satisfy the boundary condi- 
tions at the 58 corners H of the cylinder, i*e*, at r * R and z = H/2* The 
problem can be solved by numerical methods using relaxation techniques 
but is rather complicated and will not be used here* 

An approximate method of solution is to treat the reactor as having 
a radial reflector and bare ends but having a height of H + 2^ where % 

is the end reflector savings given by 

<r _ L , - 1 ( PcjjLr 1 , JjL \ 

~ 6 V Dr tcmh Lr / 

Therefore in the flux distribution expression, 



fi* - 



TT 



Since the reactor core is considered to have only a radial reflector 



where 



R+Sr 



< _ L t -l f OcBLr , , Tr 'S 
" q wn ^ — tanh ) 

From Appendix A it can be noted that the thickness of the reflector is 



uniform about the reactor core so that 



T r = T, * T 

From this, it follows that 



Following this preliminary discussion, the development of the reactivity 



expression follows closely that of the extrapolated dimensions case, 
only consideration is directed toward the effect of the reflector savings 
concept rather than the extrapolation distance. The resulting expression 
for the reactivity for this case is 



F / tP_ + 1C + l__ %£a _ L x O 1>0 \ 

Z.O V P C + <*,0 £« " KooC D J 



where 



± • 2ff(h»*h) 

H ■* Z% Z ' H + 2$ 



J_ 

Z 






stth 

H 4 



G = 



TTH i _ > _ aiTH 

• + — 5 1 n 

stCH + a'S) *i H + 



and 



= 



X<**- ( IT \ x S 

g p 



S = 



TT 

H+21 




6 t'm 



atrh 

H*a-s 



_L <; n 

z- 
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IV. COMPUTER ANALYSIS 



A. DERIVATION OF NUMERICAL COEFFICIENTS 

In the preceding sections, equations for the change in reactivity 
of the perturbed reactor system for three different models were developed. 
These equations were: 

Bare core geometry* 



Extrapolated dimensions* 



f - JL ISl i 50. \ 

J Z.A \ P + C ~ KcoC S Q “ £<*,0 D ) 



Reflected core* 



f - Jt-C<L „ H. 4 - _L_ SL') (22- 

j Z» V P C S f K„C X D ) ' • 

As can be seen, each of the equations contains the terms P, C* 2^, K 



as a whole, are constants, and are evaluated at their initial values* 
prior to the reactivity drop. From Appendix D, where expressions for 
these parameters are developed and numerically evaluated using data 
from the AGN»201, the following results are obtained. 




If. ISif L_ lis. _ If 

P C Sf <coC £(* ” <«C D 



) 



( 20 ) 




P = 0.9473 



C = 0.7414 
C s 0.06507 cm 
K K - 1.510 
£ - 0.1010 cm” 1 



a 




D * 0.2)99 err 



Before these reactivity equations describing the three models can 
be used, a discussion of the term must be presented. Recalling that 

c s 6 ^ 0r| \tJ 



where B is a constant and 



< - l — 

0 

the problem is to determine the manner in which varies as a func- 
tion of the moderator density. The Fermi Age, *C e , is proportional to 

the inverse of the product of the logarithmic energy loss per collision, 
f , and the macroscopic scattering cross-section, 

t cc — 

C# f 

To make this relation an equality, let the constant of proportionality 
be 6 t , so that 



to 









From Appendix D s 
and 

Therefore « 



Co - 8.459 cm" 
f = 0.9090 



= 1.516 



cm 



17.75 



Since <£>3 is a function of the moderator density, the relation between 



t'o and hence C, and 



is 



- 



17.15 



1 • Ba re Cylindrical Geometry 

The term e in the bare core geometry equation (20) has been de- 
fined as 



From Appendix D, 
and 



r / TT \ Z Y 

e - « ' (1 r) ‘q 

= 0.03587 cm ' 2 



H = 12.7 cm 
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Y 

Therefore O = 0.03587 + 0.1695 (23) 

Substitution of the values for the core parameters and (23) into equa- 
tion (20) yields the numerical equation for the reactivity for the bare 
core geometry model, which is 

$= -|-[(/.05C.)SP + * (15.37)^- 

- (0.3I7Z+ O.lW-^ID] < 24 > 



2. Extrapolated Dimensions 

In the equation for reactivity in the extrapolated dimensions 

model, (21), b represents the quantity 

, X«* /TT N* - W 

b “ T + VH' / T (25) 

From previous analysis it has been assumed that d, the extrapolation 
distance, can be calculated from the equation 



d = 






From Appendix D, 2^ * 1.516 cm 

so that d 53 0.47 cm 

The extrapolated dimensions are therefore 

R° - R + d » 13.17 cm 
H° - H + 2d « 25.07 cm 

In equation (25), for the extrapolated dimensions case 



ZACS .i 

56 — - 0.1826 cm 

R' 

X « J,(«xR)- J o 0*K) Jj, C«<i) 
and E = + J, Z R) 

From a table of Bessel functions, and since cX R = 2.319, 

J 6 0*R) * 0.04531 
J ( (*R) - 0.5364 



3 7 



J, (<*R) - 0.4173 



So that X = 0.5175 
and E = 0.2898 
Substituting these values into (25) yields 



0,05954 + 0.01631 - 



w 



( 26 ) 



Inserting the values for the core parameters and the expression for b, 
(26), into equation (21 ) 9 gives the numerical equation for the re- 
activity in the extrapolated dimensions model as 

S' 5c[(i.05C.)^P +<I.3M<0SC + < ,5 -37)*52 4 - ( 52*. 



- (o.5Z<>5 + O. |SS2L^-)*D] 



(27) 



3 . Reflected Core 

To evaluate the numerical coefficients in the reactivity 
equation for the reflected core model* (22)* the first step is to 
evaluate the reflector savings term* ? 

-g- ton"' C -tcrnh ) ( 28 ) 

From Appendix A* T = 20 cm 

The values for B and D have already been given as 

B = 0.2298 cm" 1 
D * 0.2199 cm 

Since the reflector material is high grade graphite, 

L = 54 . 2 cm 
r 

D = 0.94 cm 
r 

Therefore, evaluation of (26) yields 

% = 3.49 cm 
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so that 



R + §' = 16.19 cm 



H + 2 % - 31.11 cm 

For the reflected case, the parameter OC in the expression for the flux 
distribution has been given as 





2.^05 „ -1 

a — ^ - 0.1485 cm 


so that 


<X R = 1.887 


Since 


X = JoWJ^R) 


and 


E * J 0 ‘(-R)+ 


Substitution of the values 


for the Bessel functions of 



J 0 (<*£)= 0.2898 

J.faR) = 0.5814 
J^R)- 0.3269 

yields X * 0.4868 

E - 0.4218 

Substitution in the expression 

*<**■ ( 7T \ l 5 

a= E * V. H + ) F 

gives an expression for a of 

a = 0.02597 + 0.01020 -|- ( 29 ) 

Putting this expression, (29), and the values of the core parameters 
into (22) gives 

$ * TZ [(l.oeO^P + 0.3^)1C + 05,37)^ - C8.VW)*82 a 

- (o.ZZ9? + 0.090*0 oo) 

as the numerical equation for the reactivity in the reflected core model. 
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B. FLUX WEIGHTING OF CHANGE IN MODERATOR 



It will be noted that four-figure accuracy is maintained in the 
evaluation of all the above parameters. This is not strictly justified 
on the basis of the accuracy of the data available on the AGN-201 re- 
actor, but was made necessary by the very small magnitude of the reactiv 
ity change, 0,00073, which was to be accounted for. The validity of the 
above numerical results is strengthened somewhat by using a weighted dis 
tribution of the change in density within the fuel core. It seems al- 
most certain that the change in hydrogen density is not uniform within 
the fuel core, but is instead higher at the center, where the neutron 
flux is highest and decreases towards the edges of the core. While the 
density change was assumed uniform within a single disk within the core, 
the distribution of the density change in the vertical direction was 
^weighted’’ by a cosine function corresponding to the shape of the flux 
distribution in the vertical direction. 

t fa 

Define a weighting function s cos 

H 

Using the distance of the centerline of each disk (z^) from the zero 
ordinate of the z-axis of the fuel core, the fraction of the total 
density change within a single disk A n 9 can be found from the expres- 
sion 



C03 ~ 






10 

2 C,o5 

n * l M 



Z .ci n 



The results of applying this equation to the data listed in column two 
of Table 1 is shown in column four of that same table. 
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TABLE 1 



Disk 


"d 


cos 


A 


1 


4.66 


0.02973 


0.00584 


2 


4.37 


0.1253 


0.02464 


3 


4.00 


0.2455 


0.04828 


4 


3.44 


0.4198 


0.08256 


5 


2.68 


0.6323 


0.1243 


6 


1.94 


0.8012 


0.1576 


7 


0.83 


0.9626 


0.1893 


8 


-1.00 


0.9458 


0.1860 


9 


-2.50 


0.6773 


0.1332 


10 


-4.00 


0.2455 


0.04828 



C. COMPUTER PROGRAMS 

All of the discussion in this section to this point has served to 
establish the basis for the final application of perturbation theory to 
the analysis of the behavior of the AGN~201 reactor,, Appendix E con- 
tains a digital computer program, the results of which demonstrate the 
reactivity change of the reactor as a function of moderator density change 
caused by operation of the reactor* In the program, the change in moder- 
ator mass starts at 10 gm and increases to 50 gm, in increments of 10 gm. 
This overall change in mass is weighted according to the results listed 
in Table 1, to find the density change per disk* For each change in 
moderator mass, the program uses the reactivity equations (24), (27), 
and (30) to calculate the reactivity drop per disk and for the entire core 
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for each of the three models describing the reactor core. The numerical 
results of this program are presented formally in the next section. 
However^ at this points it is sufficient to state that the density 
change predicted by the perturbation theory equations corresponding to 
the actual value of the reactivity drop equal to 0.073% are as follows: 

Bare core geometry: = 17.9 gm 

Extrapolated dimensions: &Mp = 15.1 gms 
Reflected core: AMp “ 13 o 8 gm 

These values for the total core mass change are used in the computer pro- 
gram presented in Appendix F, where the effect of adding a new fuel disk 
of known properties is analyzed. Using the data for the new fuel disk 
listed in Appendix A and the results of the first computer program to 
determine the present properties of the old disks currently in the re- 
actor, the resulting excess reactivity of the reactor is predicted when 
the new disk displaces a disk of similar size by again using the numeri- 
cal perturbation equations (24 ) 9 ( 27 % and (30). The results of this 

program are presented and discussed in the following section. 
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V. RESULTS AND CONCLUSIONS 



The numerical results of Computer Program 1 presented in Appendix E 
are listed in Appendix G and plotted graphically in Figure 2. It ap- 
pears from Figure 2 that a linear relationship exists between the reactiv- 
ity of the reactor core and the change in moderator density within the 
core. Taking the slopes of these plots for each of the three models shown, 
the following expressions can be derived : 

Bare core geometry? « -0.00004095 M 

P 

Extrapolated dimensions? Al = -0.00004888 

Reflected core? Af = -0.00005302 M 

P 

Also, from Figure 2 it is possible to obtain the change in polyethylene 
mass for each model which corresponds to the observed change in reactiv- 
ity of the reactor of 0.00073. These changes are 

Bare core geometry: M - 17.9 gms 

P 

Extrapolated dimensions: = 15.1 gms 

Reflected core: M = 13.8 gms 

P 

In previous discussion a value of the weight of hydrogen bled off from 
the reactor during its lifetime of operation was estimated as 23.14 gm. 
While this value is slightly higher than that predicted by the perturba- 
tion theory equations, no true error analysis can be conducted. The 
value of 23.14 is actually only a crude approximation in itself, but the 
fact that this experimental value and the theoretical values are reason- 
ably close to one another lends validity to the results of this analysis. 

It can also be stated that the average flux level of 20 watts used to cal- 
culate the value of 23.14 gm must be close to the actual value of the 
reactor flux. Since this is the case, the arguments discounting fuel 
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Reactivity Loss (%) 



0.0030 




Decrease in Moderator Core Mass (gm) 
Fig. 2 

Reactivity Loss versus Decrease in 
Moderator Core Mass 
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burnup and reactor poisons as sources of reactivity drop are reinforced. 

Interpolation of the data in Tables 6 and 7 of Appendix G to corres- 
pond to the above changes in moderator mass allows an analysis to be 
made of the reactivity drop in the reactor core as a function of vertical 
position. These results are plotted in Figure 3, As can be seen, the 
curves for the three models are different. The extrapolated dimensions 
model curve, in general, shows a higher overall value for the reactivity 
decrease than the bare core model curve. The difference noted between 
the reflected core model curve and the extrapolated dimensions model 
curve is that the reflected core curve is slightly M f lattened® 8 . The 
reactivity decrease is higher at the top and bottom of the core and lower 
in the middle for the reflected core. The variations among the three 
curves reflect the various shapes of the assumed flux within the core 
for each model. 

Using the density changes for the three models of 17.9, 15.1, and 
13,8 gm to determine the present core density of the reactor, the re- 
sults of Computer Program 2 in Appendix F show that the effect of inter- 
changing the new disk with an old disk of the same size, (either disk 
4,5, or 6), results in a reactivity increase which exceeds the licensed 
limiting value of 0,40%. This can be confirmed by referring to Table 2 
or Appendix H. The large reactivity increase can be attributed to the 
large increase in moderator density between the old disks and the new 

3 

disk. This increase in density is approximately 0.2 gm/cm higher for 
the new disk. By modifying Computer Program 2 to print out the values 
of ^C and $D, it was found that these two positive terms in the 
reactivity equation, ( JD is negative, but is multiplied by a negative 
coefficient) which are functions of moderator density, increase substantially 
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Reactivity Decrease (7.) per Unit Disk Thickne 




Fig. 3 

Spatial Variation of Reactivity Within Core 
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REACTIVITY 



\. MODEL 
DISK 


Bare 

Core 

i 

1 


Extrapolated 

Dimensions 


Reflected 

Core 


4 . . 

(Z d =3.44) 1 

i 

\ 

l 


0.004392 


0.004816 


0.005767 


N 

CL 

II 

ho Ln 
0 

00 


0.005300 


0.005938 


0.006643 


i 

6 

(Z d =1.94) 

i 

1 

i 

i 

i 


0.006481 


0.007351 


0.007707 



TABLE 2 



when the new disk displaces an old disk. These increases more than com- 



pensate for the negative contribution to reactivity caused by the de- 
crease in fuel density of the new disk as compared to the old disks. 

For purposes of illustration, consider the values obtained from Computer 
Program 2 listed in Table 3, These values are those obtained if the new 
disk were to displace Disk 4 using the reflected core model, 

TABLE 3 
% P = 0.006269 

- 0.04459 
$<T + - 4.931 b 

=-0.00003234 cm' 1 
5.892 b 

< 8<J C * = 0.02357 b 

- 0.004563 cm" 1 
1)2* - 0.004523 cm" 1 
2D =-0.02916 cm 

From this table it can be seen that the changes in parameters caused 
directly by moderator density changes, when multiplied by their corres- 
ponding coefficients in the reactivity equation, are of about the same 
order of magnitude as those caused indirectly by the change in cross- 
section caused by the change in neutron temperature discussed earlier* 
Thus, it is concluded that a simple interchange of the new fuel disk with 
any one of three similar size disks which are presently in the reactor 
core will produce too large a value for the excess reactivity. 

The conclusions drawn from this perturbation theory analysis of 
reactivity may be itemized as follows! 
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1) The decrease in reactivity due to loss of moderator within 
the core is a linear function of the change in core density* 

2) Due to operation of the AGN-201 reactor over a period of 
several years, the moderator mass has decreased from between 14 to 18 gms. 

3) The interchange of the new fuel disk with an old disk of 
similar size currently in the reactor and in its present position within 
the core will produce an excess reactivity exceeding a value of 0.40/L 
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VI o RECOMMENDATIONS 



There are several areas into which further study of the problem of 
the reactivity analysis of the AGN-201 reactor could be directed. Based 
on the results obtained in this thesis 9 a direct displacement of an old 
disk of the same size by the new disk would exceed the limit of 0.40% ex- 
cess reactivity. The perturbation equations developed here could be used, 
however, to predict the change In reactivity resulting from any new combina- 
tion of disks in the core, such as displacing an old disk with the new 
disk but removing the smallest top disk, adding the new disk at the top 
of the core, etc. Only slight modification of the data input into the 
computer programs would be necessary to permit this. In other words, the 
equations and computer programs developed here could be used to find a 
disk arrangement such that the excess reactivity approached 0.40%. 

Secondly, a chemical analysis could be conducted both to determine 
the exact products formed after polyethylene is irradiated in the core 
and the nature of products present in the evolved gases. In this thesis 
the assumption was made that only hydrogen was evolved and that no new 
type of molecule was formed in the fuel core. Perhaps a chemical analy- 
sis of irradiated polyethylene could either substantiate this or illumin- 
ate more clearly the exact nature of this phenonemon. 

Third, the method used to relate neutron temperature to moderator 
density, equation (8), was only an empirical relation which was modified 
to suit the data in this thesis. A more exact method of handling this 
problem is to treat the problem by an energy-dependent model, with the 
scattering kernel computed by the code GAKER, which is available through 
such agencies as the Argonne National Laboratory. The use of such a code, 
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if obtained and programmed to run on the IBM-360 computer would undoubted- 
ly increase the accuracy of the results of this type of analysis* 

Finally, other analytical improvements could be made in this analy- 
sis. The results of a chemical analysis, such as previously described, 
might indicate that a two-group analysis rather than one-group as used 
here, might provide better results* Also, a more exact analytical method 
of handling the reflected core model would be desirable. 
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APPENDIX A 



DESCRIPTION OF TH E REACTOR 



A. GENERAL FEATURES i 

The AGN-201 reactor is a low power nuclear reactor designed pri- 
marily for research and educational purposes. It has a fuel core con- 
sisting of ten uranium impregnated polyethylene circular disks 10 inches 
in diameter and varying in thickness from 3/16 inches to 1-1/2 inches, 
stacked on top of one another to form a cylindrical core- The uranium 
is in the form of uranium dioxide enriched to 19-93% with U-235 and is 
uniformly dispersed in the polyethylene- Inserted in the core through 
its base are four control rods- Three are the same size; the fourth is 
smaller- Two of the larger rods serve as safety rods, the third as a 
coarse control rod, and the smaller rod functions as a fine control rod. 
Since these rods are of the same material as the core, their contributions 
to reactivity considerations are handled by incorporating their properties 
into the disk of which they form a part when fully inserted into the core. 
This is both justified and necessary when calculating the reactivity in- 
crease caused by adding a new fuel disk since insertion of the rods in- 
creases reactivity and the maximum excess reactivity value is not to be 
exceeded. 

Surrounding the fuel core on all sides is a reflector of high density 
(1-7 gms-cm ) graphite, 20 cm in thickness- Shielding is provided by a 
10 cm thick lead shield, a thick steel tank having a radius of 43.4 cm, 
a borated-water-f illed shield tank 6-1/2 feet in diameter, and finally, 
a housing of large concrete blocks which encloses the entire reactor to 
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provide protection against radiation while operating at higher power 
levels. The general features of the AGN-201 reactor are shown in Figure 
4. The core tank of the reactor is fluid tight and leak tested to a pres- 
sure of 5 psig. To handle pressure build-ups during operation and after 
shut-down, the tank is connected with copper tubing to a gas handling 
system which provides a means for measuring core pressure and a means for 
bleeding, storing, and monitoring the gases generated by the fuel disks. 

To insure against ignition of the hydrogen gas with any oxygen in the air 
which might have been initially in the fuel core, a system feeding nitro- 
gen, an inert gas, into the core tank is provided. 

B. THE FUEL CORE: 

The dimensions of the fuel core and of the ten disks composing it 
are shown schematically in Figure 5. For convenience of reference, the 
disks have been arbitrarily assigned numbers. Pertinent disk and core 
data is listed in Table 4. It should be noted that while the various 
physical properties are constant within a single disk, variations exist 
from disk to disk. 

C. NEW FUEL DISK 

Data pertaining to the new fuel disk to be inserted in the fuel core 
of the reactor is presented in Table 5. 



thermal column 




"Fig. '4 

The AGN-201 Reactor 
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1 

2 

3 

4 

5 

6 

7 

8 

9 



TABLE 4 



Fuel Core and Disk Data 



M 


V 

d 




j> P 


£235 


f 238 


292.5 


224.6 




1.196 


0.06958 


0.2796 


498.2 


482.7 




1.032 


0.05625 


0.2449 


552.4 


482.7 




1.144 


0.06648 


0.2670 


1038.1 


965.4 




1.075 


0.06259 


0.2515 


1049.1 


965.4 




1.087 


0.06313 


0.2536 


1047.2 


965.4 




1.085 


0.06303 


0.2532 


1721.2 


1930.8 




0.8914 


0.05179 


0.2080 


1722.9 


1930.8 




0.8923 


0.05190 


0.2084 


1761.8 


1930.8 




0.9125 


0.05298 


0.2129 


1784.8 


1930.8 




0.9244 


0.05371 


0.2157 




Total 


core weight s 


14821 gm 






Total 


core volume - 


12230 cm 3 






Overall core densities? 
f ^ - 0*9383 gm - cm ^ 

£ 235 = 0.05451 gm-cm" 3 

o 238 -3 

£ = 0.2191 gm-cm 





J* U = 0.2736 gm~cm 



TABLE 5 

New Fuel Disk Data 

Disk volume = 965.4 gm-cm 
Mass of polyethylene = 921.4 gm 
Density of polyethylene = 1.268 gqj-cm 
Density of U-235 = 0.06003 gm-cm^ 
Density of U-238 = 0.2287 gm-cm 
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APFEND1X B 



DISCUSSION OF THE ADJOINT OPERATOR 



The fuel core of the AGN-201 reactor is a finite cylinder* and the 
diffusion equation governing its behavior is a partial differential 
equationo While the following discussion applies equally well to such 
a case* for simplicity of presentation* consider the following differ- 
ential equation which is the one “group diffusion equation for a critical 
slab reactor shown below® 




x 




-> x 



Fig. 6 



This can be written as 



M 0 - 0 



( 31 ) 



where M is the operator 



M " dX cU + F CM 



The equation 



M + 56 « 0 



( 32 ) 



is related to equation (31) and is referred to as the adjoint equation ® 
M 4 - is known as the adjoint of M and the function is referred to as 
the importance function® In addition to satisfying equation (32)* ^ 



must satisfy the same boundary conditions as J2S * which for this case 
is to vanish at X - ~ 9 the slab half-thi ckness, The adjoint opera 

tor M4 is defined as that operator which satisfies the equality, 

r Q /t 

J-c,/ z vMw d.% = J_«y wM Net* (33) 

where v and w are any two functions which vanish at X “ “ . Sub- 
stituting the expression for M into the left-hand side of (33) and 
integrating by parts yields 

r d d n - , n | q /z- 

J-cy^, v[ 3 x DcO? + FJwdx = (vDw'-v'0^)| 

p ^ -| 

+ J.ey w bb D ax + P j V CU (34) 

v(* 70 - o 
w(i 70 - O Q 
CvDw'- v'Ow)|.^ ~ O 

f\ a /x p 

J.oy vhwdx = J.^ wMvdX (35) 

Comparison of (35) with (33) shows that for one-group theory 

M+ “ M 

When an operator and its adjoint are identicals, the operator is termed 
self-adjoint. Further* because equations (31) and (32) are homogeneous 
and Y* and satisfy the same boundary conditions* £5 and Y" 

are everywhere proportional, [5] 



But since 

and 

then* 

and from (34) 
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APPENDIX C 



EVALUATION CF INTEGRAL TERMS IN REACTIVITY EQUATION 

Evaluation of the I terns in the reactivity equation (19) is as 
follows ; 

1, = I v WVZCZf c*V 

For a cylinder*, dV r 2>TF C dr (X ~L 

and the flux distribution*, 0 , is 

0 = 00 Jo C05 



where 

and 

Therefore 



2.^08 
o< - ft 

- ft 

ir 



I,* rr-V^^VC-Oc-CaiXai.dr dz) 

Since U 9 € 9 € p and 0 O are constants within a disk and $P 9 

r 

the perturbed parameter*, is a function only of a density change which 
is assumed uniform within a single disk*, 1 can be written as 



I,® ZTT’Sf'l^C^ 0 O *T rJ 0 l (*r)c*r J h dl 



Evaluating 5 


X rJ.VOoir * TbJ/OxR) 


since 


J/C*R) • 0 




r n ^ h 4/ , 7 _ a® 

J h COS ( fit) cLl - Z.1T 


where 


^ . niii + j. . ... i . 

Q - h + z, s»n h z. 5i n 


Therefore 9 


l ( s i *i,PV<~C2f 0* R 1 HQ 


Through a similar 


analysis it is found the** 
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and 



1^ i 1>CW€P2. f 55 0 'R“ j," OR) HQ 

I, - i%Z { z^PCCR*J, l (°'R)HQ 
I H = i'»2„<Z 0 l R ! -J l l ^R)HQ 

The solution for the term is somewhat different since the diffusion 
coefficient is weighted by the gradient of the flux rather than by 
the square of the flux. The gradient of the flux can be written as 

1$ = U r br U 0 r bi 

— I c >0 v 

But, f bO ~ O 9 since it is assumed 0 ^ 0(0) 

Therefore , (V0) l = ( $ )* * ( f! T 

Taking the derivatives of 0 with respect to r and z, and squaring 
yields 

<2^ T, z (<*r) C o5 Z (^Z) +• J o 0* r ) sin 1 (fiL)] 



Therefore, X^ can be written as 



I 5 = ZR'SD2) 0 1 {” t, I t> ,i rJ, l (»>-) 0 (rJ h r ’ lS> 'oo5 l (<3Z)dZ 

_ r>ft 2 . ph't'Sh , . “) 

4 <3 J 0 r J e (<*r)dr + J h sin (02.)o(l] 



- pK _ (oCR)** 

« L rJ, C*r)dr * T" 



sR 
o 

J *h + ih 

„ cos * (<3Z)dZ = In 

I) r Jo‘ ( 0 'Odr 1 ^ 

J 



HQ 



mlh H y 

, sin^C^Octl s Tir 



where 



T TSh 
H 



V = h 4 i jin 



zr? h 



2TT (h + Sh) 



z 5in 
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Therefore 



1 5 -- W 






I 5 = { ^D0 o z R 2 J, l (*r<)H l>*Q + (-h )* y] 

The evaluation of the integral in the denomenator of the reactivity 
equation is 

V J < X'€PC2,9S 2 c(V= Znn6PC2 ( ?>/jVj 0 1 ^R)c(rJ (i M co5*(SZ)c(Z 

J* r J 0 ' (“<" ) dr = T 
J 0 Cos 1 (l 3Z) dLl - 2 

Therefore 9 

I u = \ vep 
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APPENDIX D 



NUMERICAL EVALUATION OF REACTOR CORE PARAMETERS 

A. RESONANCE ESCAPE PROBABILITY 

A general expression given by Glasstone [4] for P, the resonance 
escape probability , for a thermal -neutron reactor is 

r J_ f 2a d£ 

P(E)-- 1 exp L f J f e J (36) 



For a homogeneous reactor system, in which 2a is assumed to refer to 
the fuel only, and is taken to be independent of neutron energy in 

the resonance region 9 (36) becomes 



PCS')' «p[* 





1 

(37) 



The integral in (37) is called the effective resonance integral and 
can be represented by 




Rewriting equation (37) 



P = 




(38) 



For the AGN-201 reactor. 



v- s; * nV ♦ mV 



Dividing by N 






N 



N u N' 



<r/*< r 3 u 
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o 

However* since 0% is small and can be neglected with respect to the 
first term* 






& 



±L - JSi. 

N J 5 ' N° 



(39) 



Because 
and 

equation (39) can be expressed as 



<, e . lit? <f p 

' A p ** 



K| . £& 

N A « 



From Appendix A* 



and 



2s _ 11 (£* YV 

N° ~ f « \ A / 



J * 0*9383 gm-cm 

^ « 0*2736 gm-em 



(40) 



The evaluation of ( ~ ) is as follows* The exact mdlecular weight of 

polyethylene, A^, is unknown but can be written in terms of the (CH.) 



groups forming the molecule as 



A P * nA CH 2 - 14.0 n (41) 

Similarly* O 5 can be written as 

CRj - ncr 3 c + z.nCTj* (42) 

Dividing (42) by (41) eliminates n and hence the exact molecular weight 
need not be known co find ( )^* This same method is employed to 

find the ratio of absorption cross-section to molecular weight in a 
subsequent paragraph* Correcting the cross-sections in (42) to an 
effective neutron temperature of 340 c K through the relation 

^ <r 5 

yields 

= 3.35 b 

(T/ = /Co .“3 b 






Therefore, 



- Z..C.7 b 



Also, 



/ JiN <^ C * l*' 

' A ' A p 

A U - 270 



Substitution of these values into (40) yields 

~fjl - 2472 (43) 

An emperical relationship for the effective resonance integral in which 
U-238 is the resonance absorber is given by Lamarsh [5] as 



Substituting (43) into (44) 



- 2.73 ( ) 



x 0.486 



(44) 



I - 121.7 (45) 

The average logarithmic energy decrement per collision, f , for a 
mixture of compounds such as in the fuel core of the AGN-201 is given 
as 



I _ (Z a ?)% (g a g/ 
V ♦ £/ 



(46) 



* — -- fr - (T* = 0.0084 cm" 1 

where Gj * 13.7 b. 




- 0.036 

3’ < ° N <x(T') P = 1.508 cm 



The expression for ^ ^ is 



Using 



the value for ^ ^ is 



*<*. (f ^) c + 

f - 0,158 

f = 1.00 

$ - 0.914 
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Substituting these values into (44) yield 



f - 0o9090 ( 47 ) 

Inserting (43)* (45)* and (47) into (38) * the value of the resonance 

escape probability term is found to be 

F - 0.9473 



Bo FISSION CROSS “SECTION 

The equation used to calculate the fission cross-section* 2 j 



Nq * 



18 

From Appendix A* o ~ 0.05451 gm-cm 

The microscopic fission cross-section* corrected to a neutron tempera* 
ture of 34Q Q R is 

~ 466 b 0 



N ® 0.602 
a 



Since 

a 

235 

and A » 235 

the fission cross-section for the overall core is 



0.06507 cm 



C. INFINITE MULTIPLICATION FACTOR 

The expression for the infinite multiplacation factor* K ^ * is 

known as the four factor formula and is given as 

*?tPf 

The value of can be found from the equation 

7 -- -fr 

Since in this expression* the ^fuel' 1 ' refers to the mixture of U-235 and 

o P 

U-238* is given by 



p „ 









as^ 
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Therefore, 

Since 



and 



^7 



7 » v 



*3S 



2/ 235 > 2.47 



235 



J i36 N Q ~ aiS 



*36 ''ft 



Cl = 0.07764 



cm 



235 



where (T a 



556 b 



238 $ X?J% Ncv - *3* -l 

- — TTSV s 0*001236 cm 1 

A 



— 238 

where ^ * 2*23 b 

a 



- 2*038 

It is assumed that the fast fission factor, € , is equal to unity, 

based on the assumption that fast fission of U-235 does not contribute 
to the multiplication* From a previous paragraph 

P = 0.9473 

Finally, the equation for the thermal utilization term, f, is 



-f = 




2 : 



2 



u 

a 



♦ 



2 



7 



MS -1 

- 0.07888 cm 

2 0 P= s 0.02209 cm" 1 

/ \ P 

where /\ / = 0.0391 b. 

The resulting value for f is then 

f - 0.7810* 

Inserting these values into the four factor formula for K ^ 
infinite multiplication is found to be 

K ^ « 1.508 



, the 
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D. 



AB50RFTTGN CROSS-SECTION 



The total absorption cross-section for the fuel core is found by 
summing the cross-section of all the substances present in the core. 
Therefore, 



2 « 





£ 



P 

A 



a ~ 0.1010 cm 



E 0 DIFFUSION COEFFICIENT 

From an earlier discussion concerning extrapolation distance, the 
assumption was made that without introducing significant error to the 
calculations, D could be defined as 

D * jk 

where 2. * * 2, - A v ~ * f N, (t - ) 

since - 13.7 b, 

2 3 * 1.516 cm ^ 

and the diffusion coefficient equals 

D - 0.2199 cm 



F, DIFFUSION LENGTH 

From diffusion theory^ the quantity referred to as the diffusion 
length is defined as 

2 D 

' ~£7 

Since 2«<x and D have both been previously calculated 

2 * 2 
L **2.17/ cm 



6 7 



G 



FERMI AGE TERM 



The Fermi Age term must be evaluated last because it is solved in 
terms of all the other parameters in the criticality equation. Recall 
that the quantity C has been defined as 



C = 




(48) 



For a cylindrical reactor core, B is given by the relation 





B 2 - 


+ 


a 1 




where 


c< ^ 




= 0.1894 


-1 


R 


cm 


and 


- 


TT 

H 


= 0.1302 


-i 

cm 



Therefore, B = 0.2298 

By definition, 



2 = -I- u 



(49) 



The problem involved when dealing with an hydrogenous moderator, 
such as polyethylene, is that cannot be calculated exactly from 

the Fermi Age equation and approximate methods must be resorted to. If 
one were to write the critical reactor equation 

K»C . , 

1 1 " 1 

to include the built-in excess reactivity of 0.000354 of the core, the 
resulting "supercritical" equation would be 



KooC 



- 1.0004 



I n'B 1 

Since the only unknown in this equation is 0, the solution is 



C = 0.1704 



Fron (48) , 




= 0.1704. 



Since B is known, this is an implicit function in terms of 
iterative solution yields 




An 
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Therefore from (47) 9 



£. 0 = 0oI985 cm 
t 0 « 8.459 cm 2 . 



H. SUMMARY 

The results of the preceding paragraphs may be tabulated as 

P = 0.9473 
£ f - 0.06507 cm' 1 
K >00 =■• 1.508 

- 0.1010 cm" 1 

D = 0.2199 cm 

2 2 
L - 2.177 cm 

C “ 0.1704 

- 8.459 cm 2 



I 
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APPENDIX E 



COMPUTER PROGRAM 1 



COMPUTER NOTATION 

A5 = at neutron temperature of disk 

A8 = at neutron temperature of disk 

A235 = atomic weight of U-235 
A238 = atomic weight of U-238 

235 

AC5 - tfo at neutron temperature of 340°K 

AC8 = at neutron temperature of 340° K 

AN = Avogadro*s number 
APA = (4-)" 

AU02 = atomic weight of uranium dioxide 
B = buckling 

CCD = C for a fuel disk with change in moderator mass 
CD = C for a fuel disk 
DD = D for a fuel disk 
DEAP = 1)2/ 

DELC = 

DELD = ** D 

DELEA = 

DELEF = 

DELH = disk thickness 
DELP = 
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DELPLC = change in density of polyethylene within the fuel core 

DELPLD = change in density of polyethylene within a fuel disk 

7 j £ 

DELSA5 = ■VTc, 

< r* 

DELSA8 = 

DELSF = ^ 

DELTA = A 

DENP = density of polyethylene in a disk 

DENPLD » density of polyethylene in a disk with change in moderator mass 
DENU02 = density of uranium dioxide in a disk 
DENU5 = density of U-235 in a disk 

DENU 8 = density of U-238 in a disk 

DDD = D for a fuel disk with change in moderator mass 
EAP = 

EAU « 

EO = t. 0 

EPA = 

ERI = effective resonance integral 
ESP = £3 

EST = t s 

ESU02 = t $ 

GA3 - non-l/v absorption factor for U-235 
GA 8 = non-l/v absorption factor for U-238 
GF = non-l/v fission factor for U-235 

H * distance of nearest edge of fuel disk above or below the zero 
ordinate of the Z-axis of the cylindrical fuel core 

HH = height of fuel core 

HP = extrapolated height of fuel core 

HR - equivalent height of reflected fuel core 
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PD = P for a disk 



PM = mass of polyethylene in a fuel disk 

PPD = P for a disk with change in moderator mass 

RHOD = reactivity of fuel disk 

RH0R1 = reactivity of fuel core for bare core geometry model 

RH0R2 = reactivity of fuel core for extrapolated dimensions model 

RH0R3 = reactivity of fuel core for reflected core model 



SA50 = 


0) 

rt 


neutron 


SA80 = 


% at 


neutron 


SABAR5 


- ?.”• 




SABAR8 






SF = 






SFBAR 


II 




SFO = 


at 


neutron 


SIBAR 


II 




SIP = 


f e 




SIU02 


II 

c 




SN = 


St 





temperature of 293° K 
temperature of 293° K 



temperature of 293° K 



SPA 



P 



suo2 = <r 3 u 

TEFF = effective neutron temperature 
TM = moderator temperature 



TO = 

VD = volume of disk 
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APPENDIX F 



COMPUTER PROGRAM 2 

A. COMPUTER NOTATION 

The majority of the symbols appearing in this computer program 
are the same as in Computer Program 1 and have the same meaning. 
Therefore, reference should be made to Appendix E with the exception 
of the following terms. 



DEN PND 


= density of polyethylene in new fuel disk 


DNU5ND 


- density of U-235 in new fuel disk 


DNU8ND 


= density of U-238 in new fuel disk 


DUO 2ND 


® density of uranium dioxide in new fuel disk 


RHOMAX 


= value of excess reactivity of the core after interchanging 
old and new disk 
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APPENDIX G 



RESULTS OF COMPUTER PROGRAM 1 
TABLE 6 

Change in Moderator Core Mass = 10. gm 



Disk 


Reactivity 
(Bare Core) 


Reactivity 

(Extrapolated 

Dimensions) 


Reactivity 

(Reflected 

Core) 


i 


- 0.000000 


- 0.000000 


- 0.000000 


2 


-0.000001 


-0.000001 


-0.000003 


3 


-0.000002 


-0.000002 


-0.000004 


4 


-0.000006 


-0.000007 


-0.000012 


5 


-0.000028 


-0.000034 


-0.000045 


6 


-0.000050 


-0.000059 


-0.000068 


7 


-0.000112 


-0.000129 


-0.000137 


8 


-0.000116 


-0.000134 


-0.000140 


9 


-0.000079 


-0.000091 


-0.000097 


10 


-0.000016 


-0.000019 


-0.000025 


Reactivity 


of core (bare core 


solution) = -0 


.000410 


Reactivity 


of core (extrapolated dimensions 


solution) = -0.000476 


Reactivity 


of core (reflected 


core solution) 


= -0.000530 
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TABLE 7 





Change 


in Moderator Core Mass 


= 20 , gm 


Disk 


Reactivity 
(Bare Core) 


Reactivity 

(Extrapolated 

Dimensions) 


Reactivity 

(Reflected 

Core) 


i 


- 0.000000 


- 0.000000 


-0.000001 


2 


-0.000002 


-0.000003 


-0.000006 


3 


-0.000004 


-0.000005 


-0.000010 


4 


-0.000015 


-0.000018 


-0.000029 


5 


-0.000051 


-0.000061 


-0.000081 


6 


-0.000093 


-0.000110 


-0.000127 


7 


-0.000229 


-0.000265 


-0.000280 


8 


-0.000236 


-0.000272 


-0.000284 


9 


-0.000154 


-0.000178 


-0.000191 


10 


-0.000032 


-0.000037 


-0.000048 



Reactivity of core 
Reactivity of core 
Reactivity of core 



(bare core solutior ) = -0.000818 
(extrapolated dimensions solution) = -0.000949 
(reflected core solution) = -0.001058 
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TABLE 8 



Change in Moderator Core Mass =30. gm 



Disk 


Reactivity 
(Bare Core) 


Reactivity 

(Extrapolated 

Dimensions) 


Reactivity 

(Reflected 

Core) 


1 


- 0.000000 


- 0.000000 


-0.000001 


2 


-0.000004 


-.0000004 


-0.000009 


3 


-0.000007 


-0.000008 


-0.000016 


4 


-0.000024 


-0.000028 


-0.000047 


5 


-0.000075 


-0.000088 


-0.000118 


6 


-0.000137 


-0.000161 


-0.000187 


7 


-0.000347 


-0.000401 


-0.000424 


8 


-0.000356 


-0.000411 


-0.000428 


9 


-0.000230 


-0.000266 


-0.000285 


10 


-0.000047 


-0.000055 


-0.000071 


Reactivity 


of core (bare core 


solution) = -0 


.001227 


Reactivity 


of core (extrapolated dimensions 


solution) = -0.001423 


Reactivity 


of core (reflected 


core solution) 


= -0.001587 
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TABLE 9 



Change in Moderator Core Mass = 40. gm 



Disk 


Reactivity 
(Bare Core) 


Reactivity 

(Extrapolated 

Dimensions) 


Reactivity 

(Reflected 

Core) 


i 


- 0.000000 


- 0.000000 


-0.000002 


2 


-0.000005 


-0.000005 


-0.000012 


3 


-0.000010 


-0.000011 


-0.000022 


4 


-0.000033 


-0.000039 


-0.000065 


5 


-0.000098 


-0.000116 


-0.000155 


6 


-0.000181 


-0.000213 


-0.000246 


7 


-0.000465 


-0.000537 


-0.000568 


8 


-0.000476 


-0.000550 


-0.000574 


9 


-0.000306 


-0.000354 


-0.000379 


10 


-0.000062 


-0.000073 


-0.000095 


Reactivity 


of core (bare core 


solution) = -0 


.001638 


Reactivity 


of core (extrapolated dimensions 


solution) = -0.001899 


Reactivity 


of core (reflected 


core solution) 


= -0.002119 
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TABLE 10 



Change in Moderator Core Mass = 50. gm 



Disk 


Reactivity 
(Bare Core) 


Reactivity 

(Extrapolated 

Dimensions) 


Reactivity 

(Reflected 

Core) 


i 


-0.000001 


-0.000001 


-0.000002 


2 


-0.000006 


-0.000007 


-0.000016 


3 


-0.000013 


-0.000014 


-0.000029 


4 


-0.000042 


-0.000050 


-0.000084 


5 


-0.000121 


-0.000143 


-0.000192 


6 


-0.000225 


-0.000265 


-0.000306 


7 


-0.000583 


-0.000674 


-0.000713 


8 


-0.000597 


-0.000690 


-0.000719 


9 


-0.000382 


-0.000442 


-0.000472 


10 


-0.000077 


-0.000090 


-0.000118 


Reactivity 


of core (bare core 


solution) = -0 


.002049 


Reactivity 


of core (extrapolated dimensions 


solution) = -0.002376 


Reactivity 


of core (reflected 


core solution) 


= -0.002651 
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APPENDIX H 



RESULTS OF COMPUTER PROGRAM 2 
TABLE 11 

Reactivity Change Resulting From 
Interchanging New Disk With Disk 4 



Model 
Bare Core 

Extrapolated Dimensions 
Reflected Core 



Reactivity 

0.001582 

0.002006 

0.002957 



TABLE 12 



Reactivity Change Resulting From 
Interchanging New Disk With Disk 5 



Model 
Bare Core 
Extrapolated 
Reflected Core 



Reactivity 

0.002490 

0.003128 

0.003833 



TABLE 13 

Reactivity Change Resulting From 
Interchanging New Disk With Disk 6 



Model 
Bare Core 

Extrapolated Dimensions 
Reflected Core 



Reac tivity 

0.003671 

0.004541 

0.004897 



Maximum 

Excess 

Reactivity 

0.004392 

0.004816 

0.005767 



Maximum 

Excess 

Reactivity 

0.005300 

0.005938 

0.006643 



Maximum 

Excess 

Reactivity 

0.006481 

0.007351 

0.007707 
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